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We discuss the appropriate normalization of modes required to generate a homogeneous random 
field in an open Friedmann-Robertson-Walker universe. We consider scalar random fields and certain 
tensor random fields that can be obtained by covariantly differentiating a scalar. Modes of interest 
fall into three categories: the familiar sub-curvature modes, the more recently discussed super- 
curvature modes, and a set of discrete modes with positive eigenvalues which can be used to generate 
homogeneous tensor random fields even though the underlying scalar field is not homogeneous. A 
particular example of the last case which has been discussed in the literature is the bubble wall 
fluctuation in open inflationary universes. 
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I. INTRODUCTION 

Recently, considerable attention has been directed to 
models of structure formation in an open universe jl). 
The description of perturbations in an open universe is 
more subtle than in the spatially flat case, because one 
has to make a distinction between the concept of a func- 
tion and that of a random field. The latter can be thought 
of as an ensemble of functions with a probability assigned 
to each one. In cosmology we are interested in some finite 
region around us, possibly much bigger than the Hubble 
distance. A given perturbation about the homogeneous 
background, say the density perturbation, is some func- 
tion of space-time in this region, but the specific form of 
the function is not thought to be very important. Rather, 
we make the hypothesis that it is a typical realization of 
a random field. It is assumed that the field is homoge- 
neous with respect to the translations and rotations of 
the coordinate system. It is usually also assumed to be 
Gaussian, which means that there exists an expansion 
into mode functions with independent Gaussian proba- 
bility distributions for the coefficients. From now on we 
usually take the homogeneity and Gaussianity as read, 
referring simply to a 'random field'. 

Let us first recall the situation for a scalar field, such 
as the density perturbation. In order to generate a ran- 
dom field (and also to decouple the time dependence of 
the modes in the linear approximation) the mode func- 
tions will be eigenfunctions of the Laplacian in comoving 
coordinates. In a flat universe, the most general square- 
integrable function may be expanded in terms of the 



eigenfunctions with k 2 > (as usual — k 2 is the eigen- 
value of the Laplacian). In that case these same func- 
tions also give the most general scalar random field. In 
an open universe the situation is different. In units of 
the curvature scale, the most general square-integrable 
function may be expanded using only the k 2 > 1 modes 
(sub-curvature modes). In particular a perturbation de- 
fined in any finite region can be so expanded, which is 
all we need for any cosmological application. But now 
the situation for a random field is different; to obtain the 
most general scalar random field one needs, in addition, 
the < k 2 < 1 modes (super-curvature modes). This 
fact has been known to mathematicians for half a cen- 
tury pp], though it has only recently been brought to 
the attention of the astrophysics community Q|. One of 
the objects of the present paper is to derive, in a sim- 
ple way, the normalization of the mode functions which 
will lead to a homogeneous random field (with a natural 
choice for the variances of the coefficients) . 

In cosmology we might also be interested in vector and 
(second rank or higher) tensor functions, and the corre- 
sponding random fields. In general such objects have to 
be expanded in terms of different mode functions, which 
we shall not consider. But an important special case 
arises when the vector or tensor is the spatial deriva- 
tive of a scalar. The central purpose of this paper is to 
treat this case, and focus on a fact concerning the open 
universe which does not seem to have been discussed in 
the mathematics literature. Namely, that a homogeneous 
vector or tensor random field can be constructed by dif- 
ferentiating an inhomogeneous scalar random field. We 
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conjecture that the most general random field derived 
from a scalar can be expanded in terms of the continuum 
modes plus some new mode functions with k 2 = 1 — n 2 
(discrete modes) where n is an integer. We support our 
conjecture by displaying an explicit normalization of the 
mode functions for the cases k 2 = and k 2 — — 3 which is 
suitable for generating a tensor random field. Examples 
of both cases exist in the recent literature. A discrete 
mode with k 2 — is generated by fluctuations a massless 
scalar field in de Sitter space-time (|] and the effect on 
the microwave background of such a mode was evaluated 
in Ref. ||. A discrete k 2 = — 3 mode can be generated 
by fluctuations of the bubble wall , in single-bubble 
models of open universe inflation g! 

To be more specific, we write down the mode expansion 
for a scalar random field in terms of eigenfunctions of the 
spatial Laplacian, ~S/ 2 Zki m = —k 2 Zki m - 

/OO I 
dkJ2J2 fkimZ k i m (r,6,<j>), (1) 

1=0 m=-l 

where the fkim are members of an ensemble. The pos- 
sible values of k will be discussed in what follows. We 
will consider both a continuous spectrum of modes and 
discrete values of k. We work in a spherical coordinate 
system with line element 

ds 2 = dr 2 + sinh 2 r(d6 2 + sin 2 Od(f> 2 ) , (2) 

corresponding to the homogeneous spatial hypersurfaces 
in an open Friedmann-Robertson- Walker universe. The 
normalized eigenfunctions can be written as 

Zw m (x) = n w (r)y lTn (M), (3) 

where Y[ m (9, </>) are the usual spherical harmonics on the 
two-sphere and Hki(r) are eigenfunctions of the op- 
erator 

1 d ( o d \ 1(1 + 1) 

5- — sinh 2 r — + V 2 ' . 4 

sinh r dr \ dr J sinh r 

The normalization of the Hki{r) is the subject of the 
present paper. 

In a Gaussian random field, if one uses the complex 
form of the spherical harmonics the magnitudes of 
each coefficient have independent Gaussian probability 
distributions for m > 0. The reality condition 

fkl-m^kl = ( — l)"* (fklrn^-kl)* , (5) 

which follows from F;_ m (0,0) = (-l) m Y* m (9, (/)), fixes 
the phase of the m = modes, while the other modes 
have uniformly distributed random phases subject to the 
above equation. In what follows we will find it more 
convenient to use the real form of the Yi m (9, 0); the coef- 
ficients then have fixed phases, arg(/fe; m ) = — aig(Tlki), 
and the magnitudes of fkim are independent random vari- 
ables for all m from — I to +1. 



The variance of the distribution is defined by 

(fklmfk'l'm') = A k V(k)5(k - k')5 w 5 mm , . (6) 

where the brackets denote an ensemble average. As we 
shall discuss in the next subsection, the variance is taken 
to be independent of / and m to allow us to construct a 
homogeneous random field. Several different conventions 
exist in the literature concerning the definition of the 
power spectrum V{kY corresponding to different choices 
of the prefactor Ak [0 . In what follows we shall not need 
to define a particular separation. Since the field is Gaus- 
sian, it is completely defined by the two-point correlation 
function. For the continuous case this is defined by 

<r( Xl )/(x 2 )) = J dkAkV(k) £ Z* klm ( Xl )Z klm ( X2 ), 

lm 

(7) 

and in the discrete case the integral over k is replaced 
by a sum. Here and throughout we use the abbreviated 
notation 

oo / 
lm I— m= — l 



II. UNITARITY AND HOMOGENEITY 

The above procedure generates a scalar Gaussian ran- 
dom field. The field is said to be homogeneous if the cor- 
relation function depends only on the geodesic distance 
between the two points. By covariantly differentiating 
such a field, one can obtain a tensor random field which 
is likewise homogeneous if its correlation function is un- 
affected by a coordinate change (except of course for the 
transformation that the change induces in the compo- 
nents of the tensor). It is obvious that differentiating a 
homogeneous scalar random field always gives a homo- 
geneous tensor random field, but we shall demonstrate 
that in some circumstances one can also obtain a homo- 
geneous tensor field by differentiating an inhomogeneous 
scalar field. 

The central purpose of this paper is to ask what restric- 
tion is placed on the normalization of the modes by the 
requirement that the random field be homogeneous, and 
to discuss which values of k are compatible with this re- 
quirement. The homogeneity requirement, that the cor- 
relation function depends only on the geodesic distance 
between the points, is equivalent to the requirement that 
under a shift in the origin or orientation of the coordinate 
system the Gaussianity is preserved and the correlation 
function is unaltered. We begin by showing that for a 
scalar field this is the case if and only if such a coordi- 
nate change corresponds to a unitary transformation of 
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the basis functions.^] 

Since a coordinate change does not affect the eigenval- 
ues of the Laplacian, it will be useful to consider a field 
constructed from modes with a single value of k 



f(r, 9,(f))=y~] fklmZklm(r, 0, ' 



(9) 



lm 



The joint probability distribution for the coefficients is 



probability = N cxp — 



Him I/' 



kll. 



2A k V(k) 



(10) 



(where the normalization factor N is actually infinitesi- 
mal because the sum is infinite). 

Under a change of coordinates, the mode functions un- 
dergo a linear transformation 



Z Mm (r,8,<p) = ^2u^ ml , m ,Zki' m '(r',e',(f>') , (11) 

I'm' 

and the mode expansion in Eq. (|^) becomes 

f(r, M) = £ him Y, U Lv m >ZM>rn>{r', 8', 0') 
lm I'm' 

= Yf'ki'm-Zw m '{r',e'A'), (12) 

where 



I'm' 



fkl'r 



^ Ulrnl'r 



>fk 



klr, 



(13) 



The form of the joint probability distribution, Eq. ( |10[ ) 
with respect to the transformed coefficients is clearly un- 
altered if and only if 



/ , \f'kl'm'\ 2 



^ \fklr. 



which requires the matrix U to be unitary: 



lm 



(Ul m l"m"Y ^Iml'm 1 — fit't"$n 



(14) 



(15) 



Although we have shown this only for modes with a sin- 
gle value of fc, the corresponding expression for a contin- 
uous spectrum of modes will also be unaffected, since the 
transformation does not act on k. One can show that the 
correlation function defined by Eq. (fjj), or the spectrum 
in Eq. (Q), are also unaffected for a unitary transforma- 
tion between modes. Thus the field is homogeneous if 
and only if U is unitary. 



*The invariance of the correlation function was shown in 
Ref. As in that reference we treat infinite sums as finite, 
which should be valid if the infinite sum Eq. ([?]) is uniformly 
convergent. That has been demonstrated in Ref. M for k 2 > 
0, but we have not investigated the question for k < 0. 



The unitarity requirement allows one to change the 
normalization of the modes by an arbitrary fc-dependent 
real factor, and a completely arbitrary phase. What is 
important is the dependence of the magnitude of the nor- 
malization factor on I and m. The normalization of the 
spherical harmonics Yi m (8,(f>) ensures that they trans- 
form unitarily under an arbitrary rotation (and hence 
the distribution defined in Eq. (|6j) is isotropic), but we 
have to ensure that the complete basis functions Zki m 
transform unitarily under both rotations and shifts of 
origin. 

Note that a rotation about a fixed origin leaves k and 
I fixed but mixes the different m-multipoles, while a shift 
along the 8 = axis leaves k and m fixed but mixes differ- 
ent I multipoles. Because an arbitrary shift and rotation 
of the origin can be decomposed into a rotation, followed 
by a shift along 8 = 0, followed by another rotation, the 
unitarity of the spherical harmonics under rotations fixes 
the m-dependence of the normalization. We now seek 
the correct ^-dependence of the normalization of the ra- 
dial functions Tlki(r) which ensures homogeneity under a 
shift of the origin. 



III. HOMOGENEOUS SCALAR RANDOM 
FIELDS 

Sub-curvature modes 

For modes with fc 2 > 1, we can split IIfc;(r) as 

n H (r) = JV w n H (r), (16) 

where Nki is a normalization factor, to be determined, 
and Ylki (r) are the unnormalized functions [fj"l[ 

/ _i d \ l + 1 
n fe ;W = ^(sinhr)' — - cos(gr) , (17) 



where q 2 = k 2 — 1. For these modes, q = ±Vfc 2 — 1 is 
real and the eigenfunctions are exponentially decreasing 
beyond the curvature scale [which equals unity for the 
line element in Eq. (Q)]. It is possible to choose the nor- 
malization factor Nki to give an orthogonality relation 
between the modes (not necessarily orthonormality) of 
the form 



dVZ* klm (x)Z kn , m ,( X ) = B k S(\q\ - \q'\) 6 W 6„ 



(18) 



where Bk is a finite real function of k, independent of 
I and m. This relation ensures unitarity of the matrix 
U for a coordinate shift (J], as we show in Appendix A. 
In particular, the usual normalization factor is taken to 
be HI 
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N, 



kl 



1/2 



(19) 



s=0 



corresponding to Bk = 1 in Eq. (|l8|). 

Note that the radial functions tiki (r) and the normal- 
ization Eq. ( |l9|) are real, and one can always use the real 
form of the spherical harmonics. With that convention, 
U is real and unitarity of the transformation, given in 
Eq. (|l5|), implies that the matrix is also orthogonal; 



E 

i m 



UlmV'm" ^tml'm' ~ $Vl"$m'm" 



jk 



(20) 



k — 1 modes 

For modes corresponding exactly to the curvature scale 
q = 0, the normalization factor given by Eq. ( |l9| ) vanishes 
as Nu oc q, for all multipoles. In the absence of any 
physical reason why the amplitude of these modes should 
be zero, one can construct a homogeneous field using an 
alternative normalization factor 



N U ee lim 

q^O q 



2 1 

7T n 



(21) 



which leaves these modes finite. These modes are not 
normalizable, in the sense that the corresponding Bk = 
l/q 2 in Eq. (|l8| ) diverges as q — > 0. However, orthogo- 
nality of the transformation matrix, Eq. (^0|) , still holds 
in the analytic limit, which ensures the homogeneity of 
the field generated from the fc 2 = 1 modes. 



Super-curvature modes 



For super-curvature modes the wavenumber lies in the 
range < fc 2 
inary. 



< 1, so q = ±iyl — k 2 is purely imag- 
The mode functions can be obtained from the 
sub-curvature modes by analytic continuation, giving 



fifci(r) 



_2 (sinhr)' ' 



-l_d_ 
\ sinh r dr 



i+i 



cosh(|g|r). (22) 



These eigenfunctions extend beyond the curvature scale, 
and the resulting correlation function falls away expo- 
nentially only on scales greater than l/k Q . 

These mode functions are not normalizable, and nei- 
ther are they linearly independent of the sub-curvature 
modes in any finite region (such as our observable uni- 
verse). Nevertheless, they can be used to construct a 
more general random field than is possible from the sub- 
curvature modes alone Because there is no orthog- 
onality relation as in Eq. (|l8|), one cannot fix the nor- 
malization of these modes in the way we did for the sub- 
curvature modes. However, if we use the analytic contin- 
uation of the normalization factor for the sub-curvature 
modes 



N 



kl 



s=0 



-1/2 



(23) 



which is purely imaginary for all I, the transformation 
matrix U will also be an analytic continuation. We know 
that Eq. (15) holds in the sub-curvature regime, and in 



fact it also holds in the super-curvature regime. At first 
sight it appears that we cannot appeal to the uniqueness 
of analytic continuation to demonstrate this, because the 
left hand side of Eq. (15) is not holomorphic. However, 



if real mode functions are used the unitarity relation be- 
comes the orthogonality relation Eq. (^(]), which is pre- 
served under the analytic continuation into the super- 
curvature regime. The transformation matrices U re- 
main real in this regime because all the N^i in Eq. ( [23] ) 
are purely imaginary, and so Eq. (|2^) implies unitarity 
and homogeneity. 



IV. POSITIVE-EIGENVALUE MODES 

Modes with positive eigenvalues not only extend be- 
yond the curvature scale, but actually diverge as r — ► oo. 
However we have seen that the super-curvature modes 
with < fc 2 < 1 can form a homogeneous random scalar 
field despite being non-square integrable, so one might 
ask whether a similar analytic continuation from the nor- 
malized sub-curvature modes might also give a homoge- 
neous random field for k 2 < 0. 

For k 2 < 0, the unnormalized radial functions given 
by Eq. (|2^) remain real; however the monopole normal- 
ization factor in Eq. (|2^) is purely imaginary, while the 
dipole acquires an extra factor of i and becomes real. 
Thus the normalized mode functions no longer have a 
unique phase. This means that the transformation ma- 
trices U can no longer be purely real, unless the modes 
of differing phase do not mix with one another under a 
change of coordinates. Thus Eq. (^0|), although it still 
holds, no longer guarantees unitarity, Eq. (p^[)f] 

For instance, for —3 < k 2 < 0, the I > 1 multipoles 
are all real but the monopole is purely imaginary. One 
might hope to build a homogeneous random field from 
only the higher multipoles, but the orthogonality rela- 
tion, Eq. (J2C|), involves a sum over all the multipoles. 
Moreover a distribution with no monopole in one coor- 
dinate system will in general acquire a monopole term 
after a shift of the origin, which implies inhomogeneity. 



More generally, for 1 — (n + 1) < k < 1 



(where 



n > 1 is an integer), the normalization given by Eq. (|23| ) 



'We have not been able to prove that there is no normal- 
ization for which the transformation between modes with dif- 
ferent phases becomes unitary, only that the analytic con- 
tinuation from sub-curvature modes does not give a unitary 
transformation. 
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has a phase (I — l)7r/2 for I < n, while for alH > n it has 
the same phase, (n — 1)tt/2. A shift of the origin mixes all 
the multipoles and any attempt to construct a random 
scalar field only from modes with the same phase fails to 
respect homogeneity. 

Homogeneous tensor fields when k 2 — 1 — n 2 



diverging low multipoles given in Eqs. ( pi| ) and ( p5| ) re- 
mains finite in the limit k 2 — > 1 — n 2 even though the 
matrix elements approach zero. 

However if we act on the scalar field with an operator 
which kills the lower multipoles, we obtain a homoge- 
neous tensor random field even though the underlying 
scalar field is inhomogeneous. We now discuss two phys- 
ical cases in which this does in fact occur. 



An interesting situation arises for the discrete set of 
modes when k 2 = 1 — n 2 with n > 1 an integer. These 
modes can be obtained from the discrete closed uni- 
verse modes with k 2 = n 2 — 1 by analytically continuing 
both the radial coordinate r — ► ir and the wavenumbcr 

Note that Eqs. ( p2f ) and ( p3| ) give finite expressions 
for the first n multipoles, while for the higher multipoles 
Nkiflki(r) — 0(e), where e 2 = n 2 — l + k 2 . An alternative 
is to define 



N kl = lim eJVfc,. 

e^O e 2 ' 



(24) 
(25) 



which gives finite expressions for the I > n multipoles, 
while with this normalization the lower multipoles di- 
verge as 1/e. From Eqs. ( |l7j ) and (p5|), we have, for I > n, 



life; = — r (sinh r) 



/-Id 
V sinh r dr 



l+i 



(r sinh nr) . (26) 



Rescaling all the mode functions independently of I and 
m does not change the transformation matrix U. 

We show in Appendix B that the matrix U becomes 
block diagonal in the limit e — > 0, so that the orthogo- 
nality relation, Eq. (pp|), holds separately for transforma- 
tions between the I < n multipoles and for transforma- 
tions between the I > n multipoles. 

In a closed universe it is well known that the I < n 
multipoles form a closed unitary group under coordinate 
transformations. In an open universe, while the I < n 
modes still form a closed group, the mode functions have 
alternating phases and so the transformation matrix is 
not purely real. Thus the orthogonality relation no longer 
implies unitary, and so one cannot form a homogeneous 
random scalar field. 

By contrast, the I > n multipoles all have the same 
phase in an open universe (whereas in a closed universe 
they have alternating phases). Hence the sub-matrix of 
U connecting the higher multipoles is real, and the trans- 
formation is unitary. One might think that it is possible 
to construct a homogeneous scalar random field from the 
I > n multipoles alone. However this is not possible, 
because the lower multipoles can be regenerated by a co- 
ordinate transformation. This is despite the fact that 
U becomes block diagonal, as the contribution from the 



k — modes 

As k 2 — > 0, the normaliza tion factor for the monopole 
in Eq. (^3|) gives iV o = while the higher mul- 

tipoles diverge as Nqi oc i/k. At the same time, the 
unnormalized monopole Hoo — ► 1, while the mode func- 
tions given by Eq. (^2|) vanish for I > 1 as IIo; oc k 2 . One 
can construct a finite field from the multipoles / > 1 if 
we use the functions §§0 



N i = lim kN kl , 



U M = lim 



n. 



1:1 



k^O k 2 



(27) 
(28) 



which leaves these modes finite (although the monopole 
becomes infinite with this normalization). Such a mode 
appears if one considers the quantum fluctuations of a 
massless scalar field in de Sitter space-time using an open 
universe coordinate system Q| . The anisotropy of the mi- 
crowave background sky due to curvature perturbations 
with k 2 = in an open universe has also been discussed 
recently |6j. 

The monopole NooTloa(r)Yoo(8, (f>) is a constant when 
k 2 = 0. As a result any tensor field constructed by covari- 
ant differentiation of the scalar field will be homogeneous. 
The simplest example is the vector field 



(29) 



Note that, in the notation of Ref. |12|, V' 1 is indistinguish- 
able from an intrinsically 'vector' quantity when k 2 = 0, 
as it is solenoidal ViV = 0. 



k = —3 modes 

Another interesting case of the k 2 = 1 — n 2 modes 
discussed above is n = 2. In the limi t k 2 — > —3, the 
monopole normalization is N k o — 2iy // 2/ir and that of 
the dipole is Nki = 2^/2/ 3tt, while the normalizations 
of the higher multipoles diverge as Nki oc 1/e, where 
e 2 = 3 + k 2 . On the other hand, the radial functions 



behave as n 



k-.O 



coshr, n^i 



sinhr, and n 



h oc e 



'Note there is a typographical error in Eq. (100) of Ref. @. 



5 



for I > 2. We can use the normalization given in Eqs. ( BJ ) 
and (|25| ) to render the I > 2 multipoles finite. 

To construct a homogeneous tensor random field in this 
case we need to kill both the monopole and the dipolc. 
The simplest operator which does this is 



(30) 



where 7^ is the spatial metric. The explicit form of the 
tensor components for the metric given by Eq. (ph are 
given in Appendix C. Acting on a scalar, this gives a 
traceless and (for k 2 = —3) transverse second-rank ten- 
sor. Because of the transversality, there is no distinction 
between these 'scalar' metric perturbations and intrin- 
sically 'tensor' gravitational waves, in the notation of 
Ref. Jl2| |. A physical example of this case is the met- 
ric perturbation associated with quantum fluctuations of 
the bubble wall in open inflation models |7|,|| . 

The lowest non-vanishing modes are the quadrupole 
and octopole whose radial dependence is given from 
Eq. (p6|) as 



n fe2 
n fc3 



- ( 2 cosh r 

— ( 2 sinh r 



3 cosh r 
sinh 2 r 
5 



+ 



3r 



sinh 3 ; 



(31) 



15 



15r cosh 7 



sinh 1 



sinh 3 : 



sinh 



(32) 



Note that these higher multipoles of the scalar field di- 
verge as Tiki ~ e r as r — > 00. The action of the ten- 
sor operator renders some components of finite at 
infinity, such as T rr , T r g and T r ^ which are of order 
e _r , but the remaining components of Tij still diverge 
as e r . Nonetheless, due to the form of the metric in- 
verse [7 U = diag(l, 1/ sinh 2 r, 1/ sin 2 9 sinh 2 r)] this is 
sufficient to leave scalar invariants finite at infinity, e.g. 
e" 



-2r 



V. CONCLUSIONS 

To summarize, the normalization defined by Eq. jl9| ) 
for sub-curvature modes can be used to generate a ho- 
mogeneous field because of the orthogonality relation in 
Eq. (|l8|). We have shown that it remains valid in the 
super-curvature regime, by virtue of the fact that the ra- 
dial mode functions for a given eigenvalue fc 2 have the 
same phase. 

These normalized eigenfunctions can be written, for 
both sub- and super-curvature modes, as 



where q 2 — k 2 — 1. For sub-curvature modes, these func- 
tions are real, while for super-curvature modes they are 
purely imaginary. 

If we were to multiply the mode functions by an 
Z-dependent phase in the sub-curvature regime, they 
would still be suitably normalized there but the contin- 
uation of the phase factor to the super-curvature regime 
would in general spoil the normalization of the super- 
curvature modes. An example of this would be to re- 
place the first factor in Eq. @) by T(l + 1 + iq)/T(iq). 
Both normalizations are equivalent in the sub-curvature 
regime, but the latter is not suitable for analytic con- 
tinuation to the super-curvature regime, where it gives 
an incorrect normalization and also fails to be symmet- 
ric under q — q. Although Lyth and Woszczyna Q| 
and Hamazaki et al. jjj both quoted this latter form for 
the sub-curvature modes, they did not use it directly for 
analytic continuation and in fact both these papers ob- 
tained the satisfactory super-curvature mode normaliza- 
tion given by Eq. (|3^) above. 

For fc 2 < there is no unique phase, and therefore 
one cannot use the analytically continued mode func- 
tions to construct a homogeneous scalar random field. 
However, in the specific case k 2 = 1 — n 2 , for integer 
n, the transformation between multipoles with I > n is 
unitary, which allows a homogeneous tensor random field 
to be constructed by acting on the scalar with a covari- 
ant differential operator, provided that the operator kills 
the lower multipoles. For k 2 = we have noted that 
any differential operator does this, since the monopole 
is spatially constant. For k 2 — —3 we have seen that 
the traceless symmetric second-rank tensor does this. A 
physical example of the latter is the metric perturbation 
generated by bubble wall fluctuations in open inflation 
models. The normalized eigenfunctions for the higher 
multipoles, given in Eqs. ( p4[ ) and (|2^), can be written 
as 1 



tlni(r) 



T(l + 1 + n)T(l + 1 - n) 



2 

P~ 



1/2 



-1/2 



'(coshr) 



V sinh ? 



(34) 



We recover the flat-space limit by taking r — + while 
keeping \q\r fixed. Only the sub-curvature modes survive 
in this limit, where they tend to the usual spherical Bessel 
functions Q. The super-curvature modes, which have 
I q I < 1, are not present in this limit, and the / > n 



multipoles of the discrete modes with k 2 
discarded in this limit since n — > 00. 



1 — n 2 can be 



n fci (r) 



T(l + l + iq)T(l + l-iq) 



T(iq)T(-iq) 



1/2 



P. 



-1-1/2 



iq-1/2 



(cosh r) 



V sinh 1 



(33) 
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APPENDIX A: MODE ORTHOGONALITY AND 
TRANSFORMATION UNITARITY 

Here we prove that the orthogonality of the sub- 
curvature modes implies that the transformation induced 
by shifts in origin and orientation of coordinates is uni- 
tary. An alternative proof was given in Ref. [@. This 
means that a random field generated from them is ho- 
mogeneous. 

Substituting the transformation 

Z k i m (r, M) = E U Li>m'Zkl> m >(r',6', </>') , (Al) 



into the orthogonality relation 



dVZi lm {*)Z kn -m>{*) = B k S(\q\ - \q'\)5w 5„ 



(A2) 



gives 



dVZ* klm (x)Z m , m ,(x) 

J l"m" 

E U l'm'l"'m"' Zkl"'m"' (*') 
V"m"' 

= E (tfL"m«rtf*m'!«m«iJ fc *(|?| - \q'\) , (A3) 



I" m n 



where the last equality uses the orthogonality relation, 
Eq. (|A3), for the x' coordinates. Thus we have 



5>. 



ImV' 



U l'm'l" 



(A4) 



APPENDIX B: BLOCK DIAGONALITY OF U 

FOR k 2 = 1 - n 2 



In this Appendix we shall demonstrate that the matrix 
U becomes block diagonal in the limit k 2 — > 1 — n 2 for 
integer n > 1. We split the transformation matrix U into 
block matrix form 



U 



A B 
C D 



(Bl) 



where Ai m i> m > = Ui mVm i for 1,1' < n and so on. Then 
the orthogonality condition UU T = I gives 



A A 
AC 1 
CA 1 
CC' 1 



BB T = I . 
BD T = , 
DB T = , 

DD T = I . 



(B2) 
(B3) 
(B4) 
(B5) 



where I is the identity matrix. 

We adopt the normalization of the mode functions 
given in Eqs. ( pij ) and (|2^), which leaves the I > n multi- 
poles finite but leads to the I < n multipoles diverging as 



1/e where k = 1 - 



-e . Consider a function composed 



solely of modes with k = 1 



/fc(x) 



E 



in \ u i Y) ■ 



(B6) 



If this is to remain bounded at finite r as e — > 0, we must 
have fkim = 0(e) for I < n and fkim finite for I > n. 
Under a coordinate transformation 

A(x) - E f'wrn-Nkl^kvir^Y^ie', 0') , (B7) 



where /y/ m » is given by Eq. (jl^), and we must likewise 
have f' Wm , = 0(e) for V < n. 

Thus the sub-matrix C = 0(e). In the limit e — > 0, 
Eq. (B5) becomes DD T = I and Eq. ( |b"4| ) becomes 
DB T = 0. Because D is orthogonal this implies B = 0, 



and hence the matrix U is block diagonal. Finally, we 
note from Eq. (|b|) that AA T = I. 



APPENDIX C: k 2 



3 TENSOR MODES 



and hence the transformation is unitary. 

Note that orthogonality also gives us an expression for 
the matrix Ui ml , m ,, namely 



dVZl, Vm ,{yl)Z k i m {y.) 

= j dV'Zt Vm ,{^') E C/L"m"^"m»(x') 
I" rn" 

= Uf mVm ,B k 5{\q\-\q'\). (A5) 



In this Appendix we give the actual components of the 
(symmetric) operator Ty given in Eq. ( |30| ) for the metric 
ofEq. §, 



T 

T ftl 



T, 



4>4 



dr 2 

de 2 



3 



d 



sinhr coshr— 
\ or 

d 

- sin cos 0— 
d0 



1 



sinhr V 2 



(CI) 
(C2) 



7 



sin 2 (9 sinh r (coshr^ — ^sinhrV 2 ) , (CM) 



d_ 

Or 
0_ 
Or 
0_ 

09 



cosh r \ d 
d 



sinh r 
coshr 



sinh r 
cos \ d 



sin# ) d(f> 



(C4) 
(C5) 
(C6) 



The unnormalized k 2 = —3 monopole and dipole 
modes are given, using Eq. (p^), by 



n fco (r)F O o(0,0) 

n fe i(r)y u (e,0) 

fi M (»0*i(-i)(M) 




(C7) 
(C8) 
(C9) 
(CIO) 



It is then straightforward to verify that each component 
of Tij vanishes everywhere when applied to the monopole 
and dipole modes. Thus the action of the operator in 
Eq. (pOt) on the k 2 = —3 modes can indeed form a ho- 
mogeneous tensor field from the I > 2 multipoles of the 
scalar field whose normalized radial functions are given 
by Eqs. (|| and (||). 
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